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ABSTRACT 

Most observations of the projected half-light radii of metal-rich globular clusters in a variety of 
galaxies have shown them to be ^ 20% smaller than those of their metal-poor counterparts. We 
show using multi-mass isotropic Michie-King models that the combined effects of mass segregation 
and the dependence of main sequence lifetimes on metallicity can account for this difference, under 
the assumption that clusters with similar central potentials have the same distribution of half mass 
radii. If confirmed, this would represent a new constraint on theories of globular cluster formation 
and evolution. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: star clusters — globular clusters: 
general 



1. INTRODUCTION 

The last decade has seen rapid progress in the charac- 

I terization of globular cluster (hereafter GC) systems in 
external galaxies. An important task is to disentangle the 
properties of GCs which are universal form those that are 

[ correlated with other of their properties or those of their 
host galaxies. The Hubble Space Telescope can partially 
resolve the spatial profiles of GCs well beyond the Local 
Group, and thus study their structural parameters. The 
recovery of these parameters requires modeling of the 
point-spread function, and this has been carried out us- 

[ ing different methods and instruments by various groups, 
using a range of galaxies including spirals and ellipticals 
(Kundu & Whitmore 1998; Kundu et al. 1999; Puzia 

; et al. 1999; Larsen, Forbes & Brodie 2001; Barmby, Hol- 
land & Huchra 2002; Larsen et al. 2001; Harris et al. 
2002; Jordan et al. 2004). Most of these studies have 
revealed that metal-rich (red) GCs appear to have half- 
light radii ^ 20% smaller than their metal-poor (blue) 
counterparts. 

Larsen & Brodie (2003) have advanced the only plausi- 
ble explanation so far for this size difference. They argue 
that the observed difference can arise as a projection ef- 
fect, resulting from combination of a correlation between 
galactocentric distance and size and the differing spatial 
distributions of the GC subpopulations. Assuming that 
GCs in all galaxies follow a relation between galactocen- 
tric distance and size similar to that of Galactic GCs, 
they find that this mechanism is able produce the ob- 
served size difference, albeit with some fine tuning. 

In this Letter we propose a simple explanation for 
the observed difference. We propose that the difference 
is a consequence of mass segregation and the fact that 
lower metallicity stars have longer lifetimes for a given 
mass. Assuming that the average half mass radius does 
not depend on metallicity, we model the observed light 
profiles with Michie-King multi-mass models and stellar 
isochrones, and show that a size difference of the ob- 
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served magnitude arises naturally. 

2. MODELS 

We model GCs using Michie-King multi-mass isotropic 
models. The formalism is described in Gunn & Griffin 
(1979). In what follows, we repeat some of their expres- 
sions restricted to the isotropic case to set the notation. 

Each mass class is labeled by an integer j and is 
assumed to have a distribution function of the form 
f{E) = e^^^ — 1, where E = ^v^ + ip and V is the 
potential energy. Due to the short relaxation times at 
the core, we assume thermal equilibrium there, which 
demands that Aj — (3fhj, where /3 is a constant and fhj 
is the mean mass of the j-th class. Choosing a charac- 
teristic radius rc and velocity variance Vq, and letting 
W — —i^/vQ and ^ — r/vc, the Poisson equation reads 
V|VF = —9a, where a — p/po, p = J2j Pj *he mass 
density and po = p(0). 

If we define aj — poj/po to be the fractional density 
contribution of mass class j at the center, a model is 
completely specified by a value for W at the center, Wq, 
and the {aj}. With these quantities specified, the Pois- 
son equation is integrated until W = 0, at which point 
the tidal limit is reached. The normalized mass densi- 
ties for class j, (Jj{r), are then calculated. The models 
have mass segregation, and there is no expression re- 
lating the {aj} to the total mass in class j, which is 
the quantity we would like to specify via a global initial 
mass function CC"^); defined such that the total num- 
ber of stars in the GC with initial masses between m 
and m + dm is (^(m)d'm, and a relation between ini- 
tial and present mass, m.f{m). The problem is solved 
by iteration until the {a^} and the total mass in class 
j, Mj = Jj ^f{^)C{^) dm/ Mi), are self-consistent. 
The space densities are finally projected to obtain the 
observed mass densities. 

Determining the initial mass function in GCs is a dif- 
ficult task because of the need to take into account the 
effects of mass segregation. Using Michie-King multi- 
mass models to account for dynamical effects, Paresce 
& De Marchi (2000) find that the mass functions of a 
dozen globulars are well described by a lognormal distri- 
bution with a mean TOc — 0.33 and dispersion a — 0.34. 
We will adopt this as the mass function for GCs for 
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m < O.9Af0. At higher masses, it is found that the IMF 
is well described by a power law with an exponent close to 
Salpeter (Chabrier 2003). Thus, we adopt C(™) 
for m > O.9M0. In what follows we assume a; = 2, as this 
allows us to obtain a good match to the observed size dif- 
ference in M87. While this value is certainly consistent 
with the observations, there is otherwise no fundamental 
reason for our choice. Below we comment the effect of 
varying x. 

We construct the mass classes as follows. We take 
the minimum mass in the initial mass function to be 
O.IMq and the maximum to be 30Mq. If rata is the 
mass at the main sequence turnoff, nitrgb the mass at 
the tip of the red giant branch and N — \{mto ~ 
0.4)/0.1] = (mto - 0.4)/Am, the limits of the mass 
classes are (0.1, 0.2, 0.3, 0.4, 0.4 + Am, ...,OA+{N- 
1) Am, TOto, TOtrgb, 2, 3, 4, 5, 6, 7, 8, 30). The masses in 
each class are obtained differently depending on whether 
the value of the mass is greater than mtrgb- We ig- 
nore in what follows the evolutionary stages between the 
tip of the red giant branch and a white dwarf. This 
should not affect greatly the derived light profiles, as 
those stars will have density profiles close to that of 
stars with mto, and thus will just add light to the al- 
ready dominating component. For m < mtrgb, we have 
simply Mj oc ^. mQ{m)dm. We assume that stars with 
^ntrgb < m < 8 end their lives as white dwarfs and that 
stars with 8 < m < 30 end as neutron stars with a 
mass of 1.4Mo. For the masses mtrgb < m < 8, we 
determine the relation between the white dwarf mass 
and the star's initial mass, m^di'm'), by interpolating 
linearly the data presented in Table 3 of Weidemann 
(2000). The mass of the white dwarf bins are given then 
by Mj oc mwd{'m')Cii^)dm. Neutron stars are observed 

to be born with velocity kicks of hundreds of km s^^ (e.g. 
Lyne & Lorimer 1994), and so most of them should es- 
cape from their host GCs. Of course, the presence of 
millisecond pulsars and low mass X-ray binaries in GCs 
means that some of the neutron stars must be retained. 
Here we will assume that a fraction fns = 0.05 of neutron 
stars are retained by a typical GC. The mass in the neu- 
tron star mass class is then A/„s oc 1.4/„s /g^" C(™)'^™- 
The adopted value for f„s is consistent with the results of 
a comprehensive study of neutron star retention in GCs 
by Pfhal, Rappaport & Podsiadlowski (2002). Our re- 
sults are not very sensitive to the precise value adopted 
for fns and remain essentially unchanged when varying 
fns in the range 0.01 — 0.1. 

The final ingredient to obtain the observed light 
profiles is stellar isochrones, from which we obtain 
the mass-luminosity function L{m). We used the 
isochrones (Yi et al. 2003). Observational evidence 
points to most GCs in early-type galaxies being an old 
and coeval population (see, e.g., Jordan et al. 2002). 
We will thus assume that GCs have an age of 13 Gyr, 
the mean age inferred for the GC system of M87 using 
spectroscopic line indices (Cohen, Blakeslee & Ryzhov 
1998). We interpolated stellar populations with [Fe/H] 
in {-2, -1.75, -1.5, -1.25, -1, -0.75, -0.5, -0.25, 0}, 
an age of 13 Gyr and alpha enhancement a — 0.3. Note 
that the stellar isochrones have a minimum mass of 
O.4M0. To obtain the luminosities of lower mass stars 
we used the zero-age main sequence mass-luminosity re- 
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Fig. 1. — Projected mass and light profiles for models with 
Wo = 9 and [Fe/H]= -1.5 and [Fe/H]= -0.5, typical of metal- 
poor and metal-rich GCs respectively in early-type galaxies. The 
vertical lines indicate the half light radii r^i in units of the half 
mass radius of the respective model, showing that the metal-rich 
model has a half light radius which is ~ 14% smaller than that of 
its metal-poor counterpart. 



lation of Tout et al. (1996), ensuring continuity with the 
luminosity given in the Y^ isochrones for m = O.4M0. 
With a given isochrone in hand, we determined mto, 
mtrgb and constructed the corresponding Michie-King 
multi mass model. The projected densities of each of 
the mass classes with m < mtrgb were then multiplied 
by the mean T^-band luminosity Lvj obtained from the 
isochrones as Lyj = J. L{m)Q{m) dm/ j. C{m) dm. 

As described above, the models constructed are speci- 
fied by a value of Wq and the set of {aj}. There remain 
two arbitrary scale factors, which correspond to setting 
the scale for the spatial coordinates and the overall mass 
of the system. In order to compare the models against 
each other, we set the half mass radius Thm to the same 
physical length for all models. This assumption would 
follow if the overall structure of young GCs is determined 
by processes mostly independent of metallicity, and if 
GCs were subjected afterwards on average to the same 
dynamical effects. For each model, we then recorded the 
projected half light radius rhi in units of Thm (for each 
metallicity and Wq). 

3. RESULTS AND DISCUSSION 

In Figure Q we show mass and light profiles for two 
models with Wq = 9 and [Fe/H] equal to —1.5 and —0.5. 
The half light radii are indicated in the figure, and it can 
be seen that r^i is smaller for the metal-rich model by 
14%. This is because the mass of the most luminous stars 
becomes larger as the metallicity increases, and thus their 
density profile is more concentrated. In FigureElwe show 
rhi as a function of [Fe/H] for various values of Wq. It is 
evident from the figure that for a given value of Wq , Vhi 
gets smaller as the metallicity increases. Note that the 
half light radii are roughly half the corresponding half- 
mass radii. In order to see directly the size of the effect, 
we also plot the same curves normalized to their values 
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Fig. 2. — ( Top) Projected half light radius r^i in units of the half 
mass radius as a inunction of [Fe/H]. The different curves correspond 
to different values of the central potential Wo , ranging from Wq = 5 
(upper curve) to Wo = 13 (lower curve) in steps of AW = 1. 
{Bottom) This panel shows the same set of curves normalized to 
their value at [Fe/H]= —1.5, a typical metallicity for a metal-poor 
GC. 
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Fig. 3. — The dots are the projected half light radius for a 
sample of GCs in M87 measured using ACS images. The solid 
curve represents a robust estimate of the mean half light radius as 
a function of [Fe/H], and the dashed line is the predicted behavior 
of this quantity for a model with Wo = 9 which is normalized to 
the observed value at [Fe/H]= —1.5. 



at [Fe/H]= —1.5, typical of metal-poor GCs. A typi- 
cal metal- rich GC in early- type galaxies will have [Fe/H] 
~ —0.3. The size of the effect at that metallicity is in the 
range 5 — 25%, the exact value depending on the value 
of Wo- We will take a Wq = 9 model to be representa- 
tive of a typical GC, as the light profile of such a model 
will have log(rt/rc;) ^ 1.58, where rt is the tidal ra- 
dius and we have defined the core radius of the projected 
light profiles rd in a way akin to single mass King mod- 
els, which for a concentration c ~ 1.5 is 1.05 times the 
radius at which the luminosity is half the central value. 
Thus, this concentration measure is comparable to a typ- 
ical concentration of a Galactic GC (Harris 1996). From 
Figure[21we see that a typical metal-rich GC in an early- 
type galaxy will be observed to be 20% smaller than 
a corresponding metal-poor GC. This is consistent with 
the observations, and thus our models can explain them 
by assuming that GCs have universal physical properties 
and the combined effects of mass segregation and stellar 
evolution. 

We can go beyond the mean difference in size, as our 
procedure gives a definite prediction for the behavior of 
the half light radius as a function of [Fe/H]. The proper 
way of comparing with the observations would be to in- 
put the distribution of Wq, and then get the predicted 
behavior for rhi ■ As we don't have the distribution of cen- 
tral potentials available to us, we will assume as above 
that the population of GCs is well represented by a model 
with Wq =9. In Figure 13 we show the measured half 
light radii (average of (7475 and 2475 measurements) with 
uncertainties less than 0.5 pc for GCs in 1VI87, using data 
from the ACS Virgo Cluster Survey (Cote et al. 2004). 
These measurements form part of a systematic investiga- 
tion of structural parameters for GCs in Virgo galaxies 
(Jordan et al., in preparation). The solid line is a robust 
smoothing of the data done with the Lowess (Cleveland 



1979) method. While we show in the figure only GC can- 
didates with 1 pc < r/i; < 4 pc, there was no restriction 
in Thi for the analysis. The values of [Fe/H] are obtained 
from (3475 — 2475) as described in Jordan et al. (2004). 
The dashed line is the predicted behavior of our models 
for Wo — 9, where we have set the normalization such 
that the curves coincide at [Fe/H]= —1.5. The agreement 
is very good, especially when considering the crudeness 
of comparing with a single value of Wq. 

The results above depend on the assumption that, 
given a value of Wq , the average half mass radius does not 
depend on [Fe/H]. This assumption is appealing in that 
it points to a universality in the formation and evolution 
process of GCs. If clusters observed today with a certain 
central potential were formed with the same average half 
mass radius and were subjected, on average, to the same 
dynamical effects from the potential field of their galaxy 
and internal mass loss processes, they should have on av- 
erage the same half mass radii. Individual GCs might of 
course been subjected to quite different histories. The 
process of GC formation is yet to be fully understood, 
and thus there are few theoretical handles that would let 
us assess the plausibility of assuming a constant average 
half-mass radii with metallicity. Some proposed forma- 
tion mechanisms determine the overall scale of the proto- 
GCs by mechanisms that should be largely independent 
of the metal content, such as cosmological reionization 
compression of subgalactic halos (Cen 2001), formation 
out of dense cores of supergiant molecular clouds (Harris 
& Pudritz 1994) or on the high mass and pressure clouds 
of gas partitioned by supersonic turbulence (Elmegreen 
& Efremov 1997). Observationally, Larsen (2004) finds 
no evidence for variations on the average sizes of young 
stellar clusters in a sample of nearby spirals. Overall, 
the assumption seems certainly plausible in light of our 
current understanding of GC formation. 
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There are a number of factors that can affect the pre- 
dicted behavior of our models. An obvious one is the 
shape of the mass function, as this will change the de- 
rived {aj}. As an example, changing the power law expo- 
nent for the high mass part of our assumed mass function 
to a: = 2.35 reduces the magnitude of the difference by a 
factor of ~ 0.86. Another factor which should have an ef- 
fect insofar as it will change the evolutionary state of the 
stellar populations is age. If the populations are coeval, 
changing the age from 13 to 11 Gyr reduces the magni- 
tude of the difference by a factor of ~ 0.86. If there is 
an age difference between the subpopulations they would 
have evolved dynamically for a different total time and 
the assumption of them having the same half mass radius 
would be less warranted. Predicting the effect of varying 
the age properly would thus require knowledge of how 
the average half mass radii evolves through dynamical 
effects. This caveat notwithstanding, the fact that most 
models show rhm to be a rather stable quantity over the 
cluster's evolution (e.g., Aarseth & Heggie 1998) make it 
reasonable to use the present models and the assumption 
of constant half-mass radius to get an estimate of the ef- 
fect of an age difference. If clusters with [Fe/H]> — 1 
are 3 Gyr younger than their metal-poor counterparts, 
the magnitude of the difference would be increased by a 
factor of ~ 1.5. At any rate, most determinations of the 
relative ages of GCs in early-type galaxies, where most 
of the size difference observations have been made, are 
consistent with the GCs being roughly coeval (see, e.g., 
Jordan ct al. 2002 and references therein). The relation 
mwdirn) for white dwarf remnants also plays an impor- 
tant role in determining the magnitude of the difference 
as they contribute an appreciable fraction (~ 20%) to 
the cluster mass. 

Variations of the distributions of rhm and Wq with 
[Fe/H] or galactocentric radius can potentially be im- 
portant, but the main point to be made from our results 
is that the observations can comfortably be reproduced 
without resorting to intrinsic differences in the sizes of 
the GCs as a function of metallicity. The use of Michie- 
King multi mass models is well-suited to this task, but 
for a more precise determination of the expected size dif- 
ference, and the dependence of it with variations in the 
input ingredients, it would be very useful that this ef- 
fect be followed with N-body simulations that take into 
account the effects of stellar evolution and of the grav- 
itational potential of the host galaxy (e.g., Baumgardt 
& Makino 2003). With more detailed models in hand, 
the dependencies of the predicted half light radii could 
perhaps be used with other observed properties to simul- 
taneously constraint factors such as age and the form of 
the initial mass function. 

Our models let us predict that if the GCs are roughly 
coeval, then the size of the observed average half light 



radius difference should scale with the mean metallic- 
ity of the metal-rich GCs. This is known to correlate 
with galaxy luminosity (Brodie & Huchra 1991), so the 
average size difference should scale with it. In contrast 
with the proposal of Larsen & Brodie (2003) our models 
do not predict a change in the size difference with pro- 
jected radius. Larsen & Brodie (2003) argue that some 
pointings away from the central regions in M87 do not 
show a significant size difference between the GC sub- 
populations. Although a dilution of the size difference 
with radius could be introduced in our models by ra- 
dial variations in other factor such as the mean [Fe/H] 
of the metal-rich subpopulation, it is unlikely that the 
effect will disappear at large radii without some fine tun- 
ing. Thus, a larger number of observations of size dif- 
ference at large galactocentric radii will be very useful 
in discriminating between the models. It is possible that 
the overall effect results as a combination of projection 
and mass segregation effects, the former disappearing at 
large galactocentric radius. We stress that large samples 
of GCs are needed to investigate this issue, as GCs will 
have a distribution of intrinsic radii and we need an ac- 
curate determination of the average behavior. So while 
no size difference was reported in NGC 5128 by Harris 
et al. 2002, the low number of clusters they observed, as 
they note, precludes the drawing of strong constraints. 

We suggest in light of the models presented here that 
the size differences observed so far are consistent with 
GCs having half mass radii distributions that do not de- 
pend on metallicity. If true, GCs would present us with 
another universal property, such as the shape of their 
luminosity function (Harris 2001) or their formation effi- 
ciencies (Blakeslee et al. 1997; McLaughlin 1999), which 
can hold an important clue to their formation and sub- 
sequent evolution. And even if the effects of mass segre- 
gation do not account entirely for the observed size dif- 
ference, it is clear that its contribution must be included 
when interpreting the observations and their implications 
for GC formation and evolution. 
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